Vector chiral spin liquid phase in absence of geometrical frustration 
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Making use of detailed classical Monte Carlo simulations, we study the critical properties of a 
two dimensional planar spin model on a square lattice composed by weakly interacting helimagnetic 
chains. We find a large temperature window where the vector chirality order parameter, {Kjk)~{Sj x 
Sfe), the key quantity in multiferroic systems, takes nonzero value in absence of long-range order 
or quasi-long-range order, so that, our model is the first example where, at finite temperatures, a 
vector chiral spin liquid phase in absence of geometrical frustration is explicitly reported. We also 
show that the strength of interchain interaction is fundamental in order to obtain the vector chiral 
spin liquid phase. The relevance of our results for three-dimensional models is also discussed. 

PACS numbers: 75.30.Kz, 75.10.-b, 75.40.Mg, 77.80.-e 



Geometrical frustration and/or competition between 
interactions can lead to exotic noncollinear magnetic 
thermodynamic phases, which can be characterized by 
unusual order parameters. Particularly relevant are two 
order parameters: scalar chirality {Xjki) = {^j x S^ • 
S/) and vector chirality (or spin current) {kji^) — {Sj x 
Sfc). These two chiralities present different symmetries: 
nonzero value of {xjki) implies that the time-reversal 
symmetry is broken, while parity symmetry breaking 
comes with {kji.)^0. Both of them are relevant in 
strongly correlated electron systems: a nonzero {xjki) 
gives rise to large anomalous Hall effect [T] and leads to 
orbital electric currents in frustrated geometries [5], while 
new electromagnetic phenomena emerge in Mott insula- 
tors as a consequence of induced (xjki), generated by the 
coupling between the (/«jfc) and an external homogeneous 
magnetic field ^. On the other hand, relativistic spin- 
orbit interaction leads to a coupling between the vector 
chirality and the electric polarization [3]-[7] which play 
a fundamental role in magnetoelectric properties. This 
coupling permits also to obtain experimental informa- 
tions about the vector chirality (which is difficult to mea- 
sure owing to the absence of external physical fields that 
couple directly to Kjk)'. the chiral components in multi- 
ferroic MnW04 have been detected by neutron diffrac- 
tion using spherical polarization analysis as a function of 
temperature and of external electric field [5| . The vector 
chirality, which is the argument of this letter, always ac- 
companies helical magnetic order, and it can arise from 
spontaneous Z2 symmetry breaking in systems with com- 
petitive exchange interactions [9^ , or it can be stabilized 
by the Dzyaloshinskii-Moriya antisymmetric exchange in- 
teraction in noncentrosymmetric compounds, [SI IIOI 111] . 
However, the vector chiral symmetry can be broken also 
in a magnetically disordered state. Such phase is named 
a vector chiral spin liquid phase and has been intensively 
studied in the last years. It has been predicted to occur in 
one-dimensional (Id) frustrated quantum magnetic sys- 
tems [T2HI1]- For higher dimension d it is crucial to un- 



derstand if the vector chiral spin liquid phase is stable 
also in presence of thermal fiuctuations [13]. For d=2, 
this phase has been clearly obtained at finite tempera- 
ture T by classical Monte Carlo (MC) simulation of a 
triangular lattice of spins with bilinear and biquadratic 
interactions [15) . However, for models without geometri- 
cal frustration and d=2, 3, a clear evidence of this exotic 
phase is yet lacking, even if Onoda and Nagaosa [TH [T7] , 
investigating a Ginzburg-Landau Hamiltonian describing 
helical magnets, suggest that a vector chiral spin liquid 
phase can be stabilized even in d=3. This prediction 
was questioned by Okubo and Kawamura [15], because 
the results of their classical MC simulations do not show 
any evidence of such phase, but only a first order phase 
transition to a helimagnetic order. In this context, it is 
important to note that in the quasi- Id XY organic mag- 
net Gd(hfac)3NITEt [12] (a compound with high value 
of spin, S~7/2) a 3d vector chiral spin liquid phase has 
been experimentally observed. This result is due to the 
fact that d—1 is the lower critical dimension for an Ising 
order parameter, like (kj^), so that the chiral correlation 
length, which diverges exponentially at low T, is much 
larger than the spin correlation length, which diverges 
with a power law. Taking into account the interchain 
interaction within mean field approximation a 3d vector 
chiral spin-liquid phase results at intermediate T [20j . 
However, theoretical results obtained considering also the 
interchain fiuctuations are still lacking, and a direct nu- 
merical evidence of such a chiral phase in quasi- Id system 
will be relevant. 

In this letter we present the results obtained by em- 
ploying accurate classical MC simulation techniques to 
investigate a 2d spin system composed by weakly in- 
teracting helimagnetic chains. Despite the absence of 
geometrical frustration (the model being defined on a 
square lattice), a clear separation can be observed be- 
tween the chiral phase transition temperature T^ and the 
Kosterlitz-Thouless {KT) one Tj^t separating the quasi- 
long-range spin ordered phase from the disordered one. 
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FIG. 1: (Color online) (a) Specific heat lis. T for the model 
^ with J'^O.IJ for different L (see legend), (b) MH result 
around T^ for diflferent L. (c) Maximum of the specific heat 
vs. L (error bars lie within point size). Continuum line: fit 
function, see text. 



FIG. 2: (Color online) (a) Binder cumulant for the order pa- 
rameter K vs. T for different lattice sizes; continuous lines: 
MH interpolation, (b) Binder cumulant for the helical order 
parameter vs. temperature for different lattice sizes. Inset: 
MH interpolation around Tkt. 



We consider a simple square lattice on the (x, y)-plane 
composed of N=LxL planar spins Sij (|iS'|=l), whose 
interactions are described by the Hamiltonian: 

L L 

H = — 2_^ 2-^ 1 '^I'^^J ■ "^ij+i '^ •^2'S'ij ■ Si^j+2 

i=l 3 = 1 (1) 

j labels spins along each chain, while i is the chain label. 
Intra-chain exchange interactions are ruled by a nearest 
neighbour (NN), ferromagnetic (FM) coupling constant 
Ji and a next-nearest neighbour (NNN), antiferromag- 
netic ( AFM) coupling J2 ; interchain NN spin interactions 
are ruled by the FM coupling constant J' . If the condi- 
tion (5=1 J2|/Ji>l/4 is fulfilled, the ground state corre- 
sponds to a helical order along the chains, with a pitch 
vector g||=±cos~^ (— 1/45) . In the following we take 
S = 0.3 (i.e. Ji=l, and |J2|=0.3), while T will be given in 
units Ji . Periodic boundary conditions have been applied 
along the direction perpendicular to the chains while, due 
to the incommensurate helix modulation, free boundary 
conditions were taken along the chain direction. Config- 
uration sampling has been carried on making use of the 
usual Metropolis technique, while correlations between 
sampled configurations were mitigated by microcanonical 
over-relaxed moves [5T]. For each T at least three differ- 
ent runs have been performed, each run being composed 
of 24 X 10^ MC sweeps, with the first 4 x 10^ thermaliza- 
tion steps being discarded. Near the critical regions the 
multiple-histogram (MH) methods [H] were employed. 
Results for L = 24 — 128 are reported. 

The weakly interacting chain system has been investi- 



gated by setting J'—O.IJ: In Fig.[I|i the obtained results 
for the specific heat, Cy, vs. T are reported. We observe 
a well defined narrow and sharp peak at T ~ 0.12, which 
should be ascribed to the onset of a vector chiral spin 
liquid phase. In this T range, the scaling behavior of c„ 
for different L is reported in Fig. [lb: We immediately 
observe as increasing L the peak more and more acquires 
the typical features associated with a second order phase 
transition in the thermodynamic limit. Increasing T, a 
second broad and size-independent peak is observed at 
T ~ 0.2, which is consistent with a KT scenario. In or- 
der to estimate the critical temperatures, we employ the 

Binder's fourth cumulant U4,=l-{0'^)/3{0'^) [21], where 

-I ^" 

(where 



O will be the chirality, k 

1-1 



bii X O, 



ij + 1 



K'y.. 

K'=[L{L — 1) sing||] ), or the helical order parameter, 
^TiHM, defined as mHM=K" J dq^\S{q), where S{q) is the 
structure factor, with q=(0,q\\), and the normalization 
factor K" is the reciprocal of the structure factor inte- 
gral at T=0 [55]. The Binder cumulant for different L is 
reported in Fig. [2^ and b for the chirality and the heli- 
cal order parameters, respectively. For the chirality we 
can evaluate Tk=0. 1176(6), while for the helical order pa- 
rameter, we obtain Tr-t=0. 1095(5). The data in Fig. [2| 
allows us to assert that the two critical temperatures are 
well distinguishable with (T^ - Tkt)/Tkt - 7.4%. The 
identification of the crossing temperature in Fig. [2J3 as 
the KT transition temperature can be further validated 
by making use of the finite-size scaling (FSS) relation 
XmiL) oc L'''/"; a best fit procedure gives 7/1^=1.77(3) 
(Fig. [3k), fully consistent with the KT behaviour of a 
2d planar system. Another evidence of the presence of 
two distinct critical points comes from the scrutiny of 
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FIG. 3: (Color online) (a) Logarithm of the maximum of 
XiTiHM ^s ^ function of InL. The error bars lie within the 
symbols, (b) Vortex density p vs. T~^ for Z/=128, error bars 
he within point size. Continuum line: linear regression in 
the chiral region. Dashed line: linear regression in the KT 
regime. 



FIG. 4: (Color onhne) (a) Chiral susceptibility vs. T by MH 
interpolation for different sizes. In (b) the temperature where 
the peak of Xn is located is reported vs. L. In (c) the maxi- 
mum value of Xk is reported as a function of InL. Error bars 
fall within the symbols. 



the vortex density, p [28]. In the dilute-gas approxi- 
mation, we have p~exp(— 2/i/T), where 2/i is the en- 
ergy required to create a pair of vortices 23, , and it can 
be obtained by linear fit of —In p as a function of T"^ 
(Fig. IsId) . Three different regimes can be identified: low- 
T {T<Tkt), intermediate-T {Tkt<T<T^), and high-T 
regime {T,^<T). The linear fit in the range Tkt<T<Tk, 
(solid line) shows an activation energy of dissociated vor- 
tex pairs greater than that obtained for lower tempera- 
tures T<Tkt, where all vortex pairs are bounded, with 
a clear slope change at T^^Tkt- Finally, the creation of 
other dissociated vortex pairs appears again easier in the 
region T>Tfi, where /x strongly decreases signaling the 
onset of a complete disorder [^J ■ 

A proper characterization of the vector chiral spin liq- 
uid transition involves several aspects. A first issue, con- 
cerning the order of the transition, can be coped with by 
analyzing the equilibrium energy distribution at T—T^. 
Even at the largest simulated size of the lattice, L=128, 
no double-peaked structure was observed, so that we have 
no explicit indication in favor of classifying the chiral 
transition as a first-order one. 

The universality class pertaining to the vector chiral 
spin liquid transition has been investigated by an accu- 
rate FSS analysis. In Fig. |4^ the chiral susceptibility, 
Xk, is displayed for different values of L. From the ex- 
pected dependence of the peak position temperature on 
L, Tf^{L)=Tf^ + cL~^/'^ , we can estimate the critical ex- 
ponent v, making use of the value of T^ previously ob- 
tained from the Binder's cumulant discussed above, get- 
ting z/=l. 02(5) (Fig.jljs). Analyzing the peak values of Xk 
with the FSS relation XK.{L)oiU' / '^ , we obtain the ratio 



7/i/=1.66(7) (Fig. |4};), which implies 7=1.70(8). These 
values of 7 and v are in very fair agreement with 7=7/4 
and I'—l, i.e. the proper values of the Ising universal- 
ity class in 2d. Concerning the critical exponent a for 
Cu, which for the Ising universality class in 2d is 0, the 
c^-peak values, c™'^'^(L), vs. L are very well fitted (see 
Fig. Ic) by the FSS relation proper of the 2d Ising model 
c™="(L)=yl+Bln(L)-|-CL-i,L25^. This, allows us to con- 
clude that a=0, confirming the Ising character of the vec- 
tor chiral spin liquid transition. 

We point out that the quasi-ld nature of the model is 
fundamental in order to obtain a vector chiral spin liquid 
phase in absence of (quasi-)long-range order. This has 
been explicitly checked by MC simulations we have per- 
formed assuming J'—Ji, a model investigated by Garel 
and Doniach i26j many years ago. In Fig. [sk-d a sum- 
mary of the obtained results is reported. 

For the specific heat, c^, reported vs. T, we ob- 
serve a size-independent broad peak at T~0.75, consis- 
tent with a KT scenario, while a size-dependent, nar- 
row peak at T~0.34 is found. Using the helicity mod- 
ulus T(r) [27], one is able to evaluate Tkt taking ad- 
vantage of the universal jump: T{T)/T -^ 2/tt for 
T -> Tkt- We estimate Tkt— 0-^5. For the largest 
simulated size (L=108), Xk shows a narrow peak at the 
same temperature Tk of the narrow, size-dependent, Cy 
peak and, above all, Tk is significantly lower than Tkt- 
These observations are corroborated by the data obtained 
for a new parameter defined as M=j-'^-^^mi where 



"^.=V (^ ^U ^li^ + (r ^3=1 ^h)^ i« t^e columnar 
magnetization perpendicular to the helical displacement. 
This observable turns out to be relevant both for the 
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FIG. 5: (Color online) Observables calculated for the system 
with Ji=J' for different size L: (a) specific heat; (b) helicity 
modulus; (c) chiral susceptibility; (d) susceptibility of M (see 
text). 



indications for the behavior of their 3d counterparts. As 
we move from d=2 to d=3, the KT phase transition is 
replaced by a proper second order phase transition to 
a helimagnetic spin arrangement which implies an un- 
derlying chiral one, so that it is not surprising that in 
their MC simulation Okubo and Kawamura [TB] do not 
observe a chiral phase but only the phase transition to 
the helical order, which also entails chiral order. On the 
contrary, for a Sd collection of weakly interacting heli- 
magnetic chains, it is reasonable to hypothesize that the 
phase transition to helimagnetic order occurs at a lower 
T than the chiral one, and the vector chiral spin-liquid 
phase can manifest [20] according to the scenario recently 
emerged from experiments on the quasi- Id organic high- 
spin magnet Gd(hfac)3NITEt ^. 
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chiral phase and for the establishment of the KT phase: 
Indeed, both two- and four-points correlations contribute 
to its susceptibility xm (Fig. Isji) , which displays a first 
anomaly, which progressively stabilizes at Tkt when L 
increases, signaling the onset of a quasi-order; subse- 
quently, at lower-T, xm has a second anomaly consis- 
tent with those already displayed by c„ and Xk- So, we 
can definitely estimate T^^Q.iA. We conclude that, for 
J'=Ji, we have a clear separation between the KT be- 
haviour and chiral setup, but, at variance with the quasi- 
Id case, the onset of the chiral order is established at a 
temperature T^ lower than Tkt- 

In conclusion, we have presented the outcomes of in- 
tensive MC simulations for a 2d XY classical spin sys- 
tem, defined on a square lattice, composed by weakly 
interacting frustrated chains. We observe a clear sepa- 
ration between the vector chiral spin liquid phase and 
the quasi- long-range ordered phase, with Tk>Tkt- We 
have found that, in a system without geometrical frustra- 
tion, the chirality displays a second order phase transi- 
tion, consistent with the 2d Ising universality class. This 
result confirms the intriguing possibility of an emergent 
finite-temperature phase showing chiral long range or- 
der in the absence of the helical one as investigated by 
many authors in the multiferroic context [T^HTT] . We 
found the quasi-ld nature of the system being funda- 
mental in order to observe such an exotic phase in ab- 
sence of (quasi-)long-range order: indeed, assuming the 
same NN exchange constants in both directions (J'—Ji) 
we find that the sequence of phase transitions can be re- 
versed T^<Tkt- The opposite sequence of the two phase 
transitions for the investigated 2d models can also give 
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